Numerical study of the ordering of the ± J XY spin-glass ladder 
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The properties of the domain-wall energy and of the correlation length are studied numerically 
for the one-dimensional ± J XY spin glass on the two-leg ladder lattice, focusing on both the spin 
and the chirality degrees of freedom. Analytic results obtained by Ney-Nifle et al for the same 
model were confirmed for asymptotically large lattices, while the approach to the asymptotic limit 
is slow and sometimes even non-monotonic. Attention is called to the occurrence of the SO(2)-Z2 
decoupling and its masking in spin correlations, the latter reflecting the inequality between the 
5*0(2) and Z2 exponents. Discussion is given concerning the behaviors of the higher-dimensional 
models. 



§1 Introduction 

The domain-wall method, or the stiffness method, is 
widely used in studying the ordering properties of vari- 
ous spin systems including the spin glass (SG). In this 
method, one computes by some numerical means the 
change of the ground-state energy of finite systems of 
the linear dimension L under the appropriate change 
of boundary conditions (BCs) imposed on the system 
P, 0, 13- This energy is called a stiffness energy (or a 
domain-wall energy), AE(L), which gives a measure of 
an energy scale of low-energy excitations of size L. For 
large L 7 AE(L) is expected to behave as a power-law, 
AE(L) w L , 9 being a stiffness exponent. If < 0, the 
system remains in the disordered state at any nonzero 
temperature, whereas if 9 > the system possesses a fi- 
nite long-range order at low enough temperatures with 
T c > 0. 

For complex systems like SGs in which the nature of 
the ordering is highly nontrivial, some fundamental ques- 
tions still remain in this method. First question concerns 
with the choice of the set of BCs. In principle, there 
could be various choices of the set of BCs, and the behav- 
ior of AE{L) may depend on these choices, particularly 
for small sizes practically accessible in numerical simula- 
tions. It is not generally clear which set of BCs should 
be best chosen, especially when the results depend on 
the BCs. Second question concerns with the meaning of 
the stiffness exponent. When the domain- wall energy de- 
creases with L, i.e., the stiffness exponent is negative, it 
is a common practice to relate the inverse of the stiffness 
exponent with the correlation-length exponent v associ- 
ated with the T = transition, i.e., v — l/\9\ = l/y 
ill = 1^1)0 El E| By contrast, there were reports which 
cast doubt on the validity of this simple relation . Fur- 
thermore, if the model could exhibit more than one stiff- 
ness exponents depending on the choice of the set of BCs, 



a question immediately arises which stiffness exponent 
should be chosen to estimate the correlation-length ex- 
ponent. Since, in the standard continuous (second-order) 
phase transition, only one diverging length scale is ex- 
pected, the existence of more than one stiffness exponent 
poses some problem in its interpretation. 

In fact, under certain circumstances, there could be 
more than one distinct diverging length scales, or more 
than one distinct correlation-length exponents, at a single 
continuous transition. An example of this may be seen in 
chiral transitions possibly realized in certain frustrated 
vector spin systems including the XY SG @. Frustrated 
vector spin systems often possesses a chirality degree of 
freedom due to the canted spin structure, according as 
the noncollinear (or noncoplanar) structure induced spin 
frustration is either right- or lef-handed. Chirality is a 
pseudo-scalar variable, being invariant under global spin 
rotations [50(2)] but changing sign under global spin re- 
flections [Z2) ■ Recent studies have suggested that some of 
such chiral spin systems might possibly exhibit a "spin- 
chirality decoupling" phenomenon, where the Z2 chirality 
exhibits an ordering behavior entirely different from the 
50(2) spinQ, though there still remains some contro- 
versy concerningwhether such a spin-chirality decoupling 
really occurs [EH II El 

In one possible realization of the spin-chirality decou- 
pling, the chirality and the spin exhibit two separate 
transitions at mutually distinct temperatures, whereas, 
in other possible realization, the chirality and the spin 
order at the same temperature but with the mutually 
different spin and chirality correlation-length exponents, 
u s and v K . Examples of the first class might be the or- 
derings of the regularly-frustrated two-dimensional (2D) 
XY model lal and those of the three-dimensional (3D) 
XY SG [ISllIIII- Examples of the second class the 
ordering of the regularly frustrated one-dimensional (ID) 
XY model Q and those of the 2D XY SG [HEi- A 
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firmly established example is the case of the regularly- 
frustrated ID XY model, where it has been shown rigor- 
ously that the spin and the chirality order at T — where 
the chiral correlation length diverges exponentially with 
v K = oo but the spin correlation length as a power-law 
with v s = 1 |l3j|. 

More controversial is the nature of the ordering of 
the XY SG. Some time ago, Kawamura and Tanemura 
made a numerical domain-wall study of the XY SG in 
2D and 3D 0, 0] . These authors introduced various 
types of BCs to probe the spin and the chirality order- 
ings of the model, including the periodic (P), antiperi- 
odic (AP) and reflecting (R) BCs. In particular, the 
domain-wall energy obtained under the combination of 
the P and AP BCs (P/AP), AE PtAP (L), and the one 
obtained under the combination of the R and P BCs, 
AE C (L), apparently yielded mutually different stiffness 
exponents, which were interpreted as associated with 
the spin and the chiral correlation-length exponents, re- 
spectively. These authors observed that, in 2D, both 
A£p,ap(£) and AE C (L) decreased with L, characterized, 
respectively, by mutually different stiffness exponents, 
y K ~ 0.5 and y s ~ 1.0. This observation was interpreted 
as indicating that the chiral correlation length outgrows 
the spin correlation length at the T = transition of the 
model, i.e., v K > v s 10]. These results were corroborated 
by several Monte Carlo CMC) simulations on the 2D XY 
SG [3 El El El El 113 B y contrast, on the basis 
of their domain-wall energy calculation, Kosterlitz and 
Akino claimed that the spin and the chiral correlation- 
length exponents were common at the T — transition 



In 3D, Kawamura and Tanemura observed that 
AEp,ap{L) (AE C (L)) decreased (increased) with L, 
which was interpreted as indicating that the chiral-glass 
transition occurred at a nonzero temperature, Tcg > 
0, while the standard SG transition occurred only at 
Tsg = [1(1 Ell. MC results supporting such a view were 
also reported |lllE^|. Meanwhile, later domain- wall en- 
ergy calculation by Macourt and Grempel suggested that 
AEp t Ap(L) might eventually be iterated toward strong 
coupling for larger L and that T C g > T SG > [20|. By 
contrast, Lee and Young claimed on the basis of their 
MC simulations that the spin and the chirality ordered 
at the same finite temperature Tcg = Tsg > 0, with a 
common correlation-length exponent [i| v s — v K . Thus, 
the situation remains quite controversial. 

Inspired by the numerical work of Ref.[lfJ on the 2D 
and 3D VFSG, Ney-Nifle, Hilhorst and Moore performed 
an analytic study of the ID XY SG ladder with the bond- 
random ±J (or binary) interaction |21j. The Villain's 
action was assumed there. Via the dual transformation, 
the model was mapped onto the ID charge Hamiltonian. 
Since the mapped model was still not amenable to the 
exact treatment, Ney-Nifle et al made further simplifica- 
tions, and eventually derived several analytic results con- 



cerning the domain- wall energies and correlation lengths. 
They observed that, depending on the type of the applied 
BCs and on whether the total number of frustrated pla- 
quettes is either even or odd, the domain-wall energy 
exhibits different behaviors. When the sample average 
is taken over all samples, AEj^ap is characterized by 
the chiral stiffness exponent y K — 1.899 • • •, while AE C is 
characterized by the spin-wave (SW) stiffness exponent 
y' s = 1. This is in contrast to the assignment made in 
Ref.E3 for the 2D and 3D XY SGs. 

If one looks at the spin and chiral correlations of this 
ID ladder model, one sees that the model exhibits the 
spin-chirality decoupling in the sense that there exist two 
distinct diverging lengths at the T = transition (though 
Ney-Nifle et al apparently stated otherwise), the one as- 
sociated with the Z2 chirality and the other associated 
with the SO(2) SW. The chiral correlation length £ K is 
characterized by the exponent v K = l/y K = 0.5263 
while the SW correlation length is characterized by 
the SW exponent v' s = l/y' s = 1. The full spin corre- 
lation function is the product of the Zi part with the 
correlation length £ K w T~ 526 and the SO (2) part with 
the correlation length sa T _1 . Reflecting the fact that 
the Zi chiral correlation-length exponent happens to be 
smaller than the SO(2) SW correlation-length exponent, 
i.e., v' s > v K , the full spin correlation function is domi- 
nated by the chiral exponent v K . 

The analytic result of Ref . [2]| . though quite plausible, 
is not completely rigorous. Furthermore, some of the 
results were obtained for asymptotically large lattice. In 
the SG problem, it is sometimes important, and often 
not a trivial matter, to elucidate the finite-size effect, 
e.g., how large the system must be for exhibiting the 
asymptotic large-lattice behavior. 

Thus, we feel it would be useful to perform numer- 
ical study of the ± J ID XY SG ladder in comparison 
with the analytic work of Ref.[2l]]. In the present paper, 
we undertake such numerical analysis of the ± J XY SG 
model on two-leg ladder lattices. The aim of our calcu- 
lation is threefold, (i) We wish to test the validity of 
the simplifications made in the analytic work of Ref. |2J| . 
(ii) We wish to elucidate the nature of the finite-size ef- 
fect in this ID model, (iii) We wish to further exam- 
ine the relation between the stiffness exponents and the 
correlation-length exponents in this ID model. 

The following part of the paper is organized as follows. 
In §2, we introduce the model and summarize the analytic 
results of Ref . |2*ll| . In §3, we present our numerical results 
of the domain- wall energies. The results are compared 
with those of the analytic work of Ref.j^]. Finite-size 
effects are analyzed carefully. In §4, we present our nu- 
merical results of the spin and chiral correlation lengths, 
in comparison with the corresponding analytic results of 
Ref.[21|. Relation with the stiffness exponents and the 
correlation-length exponents are examined. Finally, §5 is 
devoted summary and discussion. 
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]2. The model and some analytic formula 



the form, 



The model we consider is, firstly, the standard W-SG 
model on the ID two-leg ladder lattice with the binary 
(or ± J) interaction, whose Hamiltonian is given by 



n 



j^Si ■ Sj 

<ij> 



<ij> 



Jij coa(Bi -6j), (1) 



where Si = {S?,S?) = (cos 0*, sin 0*) (0 < 9 t < 2tt) is 
the two-component spin variable at the site i, and the 
summation is taken over all nearest-neighbor pairs on 
the ladder lattice. The site index i may be written as 
i — {x,y) with 1 < x < L and 1 < y < 2, where y = 1 
and 2 refer to the first and the second row of the ladder. 
Jij represents the random variable taking either +1 or 
— 1 with equal probability independently at each bond. 
The absolute value of the exchange interaction has been 
taken to be a unit of energy ( J = 1). 

In the following, we impose several types of BCs on the 
W-spin variables at the boundary, i.e., the periodic (P), 
antiperiodic (AP) and reflecting (R) BCs. In these P, AP 
and R BCs, we impose the relations, SiL+i,y) — S(i,y)> 



), re- 



spectively. In the R BC, we reflect the spin at the bound- 
ary with respect to the x-axis in spin space. 

The local chirality variable at the plaquette x, consist- 
ing of four spins at the sites (x. l),(x + 1, 1), (x + 1, 2) 
and (x,2), are defined by 



— y= J2 sgn(J tj )sm{ 

<ij> 



(2) 



where the summation is taken over four bonds connect- 
ing the above four sites forming the plaquette. In the 
ground state of an isolated frustrated plaquette, it takes 
a value either +1 or — 1, while in the ground state of an 
isolated unfrustrated plaquette, it takes a value equal to 
zero. Thus, the states with k = ±1 represent the two 
chiral states, with right-handed and left-handed spin cir- 
culation around the plaquette. 

Other model we consider is the effective charge Hamil- 
tonian on the dual lattice. The simplest version is the so- 
called Villain's Hamiltonian, which contains only the 2- 
body charge interaction. Although it has commonly been 
believed that the Villain's Hamiltonian becomes equiva- 
lent to the cosine Hamiltonian in the low-temperature 
limit T — ► [2l| , we have found that this is actually not 
the case: The original XY Hamiltonian mapped to the 
charge representation contains the higher-body interac- 
tions, in addition to the 2-body interaction, even in the 
T — > limit. Thus, in the present work, we also consider 
these higher-body correction terms to the standard Vil- 
lain's 2-body approximation |2fij . The explicit forms are 
given in the appendix. 

In the cases of the P and AP BCs, the two-body 
charge Hamiltonian, or the Villain's Hamiltonian, takes 



H p = £ U v m * m 3 + ^r( 2n ~ X] m * " V fi ( 3 ) 



i>3 



n 2 

Hap = Uijrriimj + —(2n - ^ m, - V + l) 2 , (4) 

respectively, where the charge variable m,-, sitting at the 
plaquette i, takes integer values 0,±1,±2... on unfrus- 
trated plaquettes, and half-integer values ±i,±|... on 
frustrated plaquettes. The variable n takes integer values 
0, ±1,±2..., while V is the "parity" variable being equal 
to zero or unity depending on whether the total number 
of antifcrromagnctic bonds on the first row (y = 1) of the 
ladder is either even or odd. The interaction between the 
charge variables Uij located at the plaquettes i and j is 
defined by 



_ ^2 e ik(i-j) 

ij ~ T ^ 2 - cos k ' 

fe 



(5) 



where the summation over the wavevector k is taken over 



±±Z. . 

' L ' L ' 



In the L — » oo limit, [/,-,■ reduces to 



Un = £=(2-^/3)'*- 



(6) 



which decays exponentially with distance \i — j\. The 
first term of Eq.(2) and (4) represents the charge-charge 
interaction, while the second term of Eq. (3) and (4) rep- 
resents the SW term, which is related to the charge part 
via the total charge JV m,-. 

In the case of the R BC, by contrast, the corresponding 
two-body charge Hamiltonian is given by 



Hn = X Uijmimj, 



(7) 



where Uij is still given by Eq.(5), but the summation 
over the wavevector k is now taken over k = ±|j, ±^S, ■ • • 
which yields Ui+L,j = ~Uij. Note that there is no second 
term (SW term) in Hr. 

To proceed furtherjNey-Nifle el al made the follow- 
ing two assumptions [2lJ. First, the charge variable m; 
is restricted to ±i on frustrated plaquettes and on 
unfrustrated plaquettes. If one labels the frustrated pla- 
quettes as I = 1,2,... Nf T , where 7Vf r is the total number 
of frustrated plaquettes, the chiral part of the Hamilto- 
nian reduces to the ID Ising Hamiltonian with Vf r Ising 
variables 07 = ±1. Second, Ney-Nifle made a further 
simplification that the charge-charge interaction, which 
originally work between arbitrary pairs of frustrated pla- 
quettes, is restricted only to the nearest-neighbor pairs 
of frustrated plaquettes. After these two simplifications, 
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the model reduces to the ID Ising chain with the random 
antiferromagnetic nearest-neighbor interaction, 

Rising = ^ Viami+i. (8) 
1=1 

The random nearest-neighbor interaction Vj > obeys 
the distribution given by 

P(V) = cV~ 1+v «, v K = 0.5263 (9) 

for smaller V , where c is a normalization constant. The 
latter simply follows from Eq.® and the fact that the 
probability to have a sequence of / successive unfrustrated 
plaquettes is given by 1/2 1 . Thus, it should be remarked 
that the distribution of the effective interactions P(V) is 
not a smooth function but a collection of delta functions 
so that Eq.(PJ must be taken with care. 

These simplifications enabled Ney-Nifle et al to spec- 
ify the explicit charge (chirality) pattern and the 
existence/non-existence of the SW excitation in the 
ground state of large enough lattices under the given BC, 
leading to various predictions on the domain-wall ener- 
gies and the correlation lengths. 

§3. Numerical results on the domain-wall energies 

In this section, we numerically calculate the following 
two types of domain-wall energies for the ID XYSG lad- 
der, i.e., (i) the root-mean square of the energy difference 
between under the P and AP BCs, AE^ap, and (ii) the 
absolute value of the energy difference between under the 
R and min(P,AP) BCs, AE C , where min(P,AP) refers to 
either P or AP BC which has the lower energy than the 
other. The following three levels of numerical calcula- 
tions are made. 

Method (A): The first method is the direct numerical es- 
timate of the ground-state energy of the cosine Hamil- 
tonian of finite L under the given BC. In calculating 
the ground-state energy, we employ the spin-quench algo- 
rithm without any further approximation, i.e., by start- 
ing from the randomly generated spin initial conditions, 
we quench the system to reach one of the local energy 
minima. These quench procedures are repeated many 
times, typically 5,000 times, until one is sure that the 
true ground state has been reached. If one goes to larger 
L, the number of local minima increases rapidly which 
makes the search of the true ground state increasing dif- 
ficult. This difficulty limits the tractable maximum lat- 
tice size to L < 35. The sample average is taken over 
10,000 (L = 35) - 200,000 (L < 15) independent bond 
realizations. 

Methods (B1)-(B3): In this group of methods, we es- 
timate the ground-state energy of the effective charge 
Hamiltonian of finite L under the given BC within the 
first approximation of Ref . [2l| . Namely, we use the Ising 



Hamiltonian but with the distant-neighbor interaction in 
estimating the domain-wall energy. In identifying the 
chirality pattern of the ground state (or the candidate 
of the ground state), we have made certain plausible as- 
sumptions, the detail of which will be given below for 
each different case: It is similar to the procedure em- 
ployed in Ref. [23, but it is systematically improved in 
the estimation of the domain wall energies by taking into 
account the distant-neighbour two-body interactions and 
also the higher-body intereactions. By the methods (Bl)- 
(B3), we examined lattices considerably larger than the 
case (A), up to L = 960. The sample average is taken 
over 100,000 independent bond realizations. 
Method (C): In this third method, we estimate the 
ground-state energy of the Villain's Hamiltonian of finite 
L under the given BC by assuming both the first and 
the second approximations of Ref.[2l|. Namely, we use 
the Ising Hamiltonian with the nearest-neighbor interac- 
tion in estimating the domain-wall energy. In this case, 
we can deal with lattices still larger than the cases (Bl)- 
(B3), and the results for asymptotically large L should 
reduce to the analytic results of Ref.|2lJ. The sample 
average is taken over 100,000 independent bond realiza- 
tions. 

We compare the results of these five levels of calcu- 
lations (A), (B1)-(B3) and (C) to examine the validity 
of the approximations made in Ref. [2l| , and to elucidate 
the nature of the finite-size effects in this model. 

In the original cosine model (1), the local chirality at 
each plaquette is given by Eq.(2). In the ground-state, 
the local chirality distribution is expected to be peaked 
around k = ±1 for frustrated plaquettes, and around 
k = for unfrustrated plaquettes. We have checked nu- 
merically that this is indeed the case. As shown in Fig.l, 
the local chirality on frustrated plaquettes takes the val- 
ues around n ~ ±1 with equal probability, while that on 
unfrustrated plaquettes takes the values around n = 0. 
Such a distribution of the local chirality enables us to 
label the chirality pattern uniquely only by the combina- 
tion of + and — on frustrated plaquettes. 

In the present ID model, as was first suggested by 
Ney-Nifle et al [2l| , the behavior of the domain- wall en- 
ergy largely depends on whether the total number of frus- 
trated plaquettes N{ x is either even or odd. We call these 
samples "even" and "odd" samples, respectively. In what 
follows, we show the results for the even and odd cases 
separately. 

§3.1 Even samples 

3.1.1 The domain-wall e nerg y: AEp ap 

For even samples, Ref.|2l| predicted that the domain- 
wall energy AE^ap was dominated for large enough lat- 
tices by the contribution of a pair of chiral domain-wall 
excitations not accompanying the SW excitation. The 
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FIG. 1: The distribution function of the local chirality k 
at each plaquette in the ground state of the ±J XY ladder. 
The lattice size is size L = 35. The data for frustrated and 
unfrustrated plaquettes are given in blue (at the center k — 0) 
and in red (near the edges k = ±1), respectively. 



chiral domain- wall may be defined here as the place where 
the "chiral overlap" between the two chirality configu- 
rations 1 and 2 under the two BCs, Oi = mi(l)m.;(2), 
changes the sign. According to Ref.jU, under the 
min(P,AP) BC, the sign of the chirality pattern alter- 
nates on frustrated plaquettes without misfit, while un- 
der the max(P,AP) BC, a pair of misfits is introduced into 
the alternating chirality pattern, but not accompanying 
the SW excitation: A pair of chirality misfits is intro- 
duced into the sample in such a way that one is at the 
weakest connection, i.e., the place where the neighbor- 
ing frustrated plaquettes are most far apart in distance, 
and the other is at the next-weakest connection satisfy- 
ing the condition m » = ^l- The latter condition is 
required to suppress the SW term in eqs.(3) or (4). The 
other possible candidate of the ground state under the 
max(P,AP) BC might be the 7T-SW state with a nonzero 
second term but without any misfit in the alternating 
chirality pattern. However, if the assumptions made in 
Ref.[21j are to be justified, the chiral domain- wall state 
always has the lower energy than the SW state, at least 
in sufficiently large lattices. 

In our methods (B1)-(B3) above, we search for the 
positions of a pair of chiral domain walls in the ground 
state under the max(P,AP) BC according to the follow- 
ing procedure. First, to specify the position of one of 
the two chiral domain walls, we apply the R BC to the 
same sample. The application of the R BC is expected 
to yield a single chiral domain wall in the sample: See 
below. By calculating and comparing the energies corre- 
sponding to all possible N& positions of a chiral domain 
wall, we determine the position of the chiral domain wall 
in the ground state under the R BC. This position is as- 



sumed to be common with the position of one of the two 
chiral domain walls introduced under the max(P,AP) BC. 
We then determine the position of the second chiral do- 
main wall in the ground state under the max(P,AP) BC 
by calculating and comparing the energies corresponding 
to all possible positions of the second chiral domain wall 
under the constraint that there are odd number of frus- 
trated plaquettes between the two chiral domain walls. 
In calculating the energy, we employ systematically im- 
porved methods (B1)(B2) and (B3) unlike the case of 
our method (C) which takes into account the nearest- 
neighbour intereaction only. In the method (BI), we use 
the two-body approximation but sum over all distant- 
neighbor interactions. In the methods (B2) and (B3), we 
take into account the higher-body correction terms up 
to the 4-body and 6-body interactions, respectively. The 
explicit forms of the higher-body terms are given in the 
appendix. 

In our direct method (the method (A) above) for fi- 
nite L < 35 samples, we have observed that, under the 
min(P AP) BC, the rule of the ground-state configuration 
of |2lJ is always satisfied, while, under the max(P,AP) 
BC, some samples obey the rule of Ref [2l|, but some 
other samples do not. In the latter class of samples, the 
ground state under the max(P,AP) BC turns out to be 
the 7T-SW state rather than the chiral domain-wall state, 
i.e., the chirality pattern completely alternate without 
misfit while the SW of a turn angle it (7T-SW) is gener- 
ated between under the P and AP BCs. An example of 
such a 7T-SW sample is shown in Fig. 2, where the spin 
configurations under the P and AP BCs are shown in the 
upper panel of Fig. 2, while the relative deviation angle 
between the spin directions under the P and AP BCs is 
illustrated in the lower panel by arrows. The appearance 
of the 7T-SW is clearly visible here. Hence, at least in a 
subset of samples of finite L < 35, the rule of Ref . [2lj| is 
violated. 

The domain- wall energy AEp^p calculated in our di- 
rect method (A) is shown in Fig. 3(a) on a log-log plot. 
The slope is estimated to be about 1.39 in the range 
L < 35, which is considerably smaller than the predicted 
value of Ref. [13], 1.899- • •, presumably due to the exis- 
tence of the SW samples characterized by the SW stiff- 
ness exponent y' g = 1. 

Then, the next question is how the rate of the 7T-SW 
samples, which breaks the rule of Ref. [2]|, varies with 
increasing L. This rate of the SW samples rsw calcu- 
lated in the direct method (A) is shown in Fig. 4(a) in 
the range L < 35. As can be seen from figure, while 
the rate decreases with increasing L for smaller lattices 
of L IS 20, it tends to increase again for larger lattices 
up to L — 35. Unfortunately, the direct calculation is 
limited to L = 35. 

The rate of the SW samples rsw is also estimated by 
using the Ising approximation (the method (B1)-(B3)) 
and the nearest-neighbor approximation (the method 
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FIG. 2: A typical example of the ground-state spin config- 
uration of even samples under the periodic and antiperiodic 
boundary conditions, between which a spin-wave of a turn 
angle -k is generated. The lattice size is L = 25. In the lower 
panel, the relative deviation angle between the spin directions 
under the periodic and antiperiodic boundary conditions is il- 
lustrated by arrows. 



(C)), and the results are also shown in Fig. 4(a) in the 
same range of L < 35. One sees from the figure that 
the rate calculated by both approximations exhibits the 
non-monotonic behavior qualitatively similar to the one 
observed by the direct method (A). In particular, the 
Ising approximation (B1)-(B3) gives the results in quan- 
titative agreement with those of the the direct method 
(A). The quantitative agreement becomes systematically 
better by including the distant-neighbour two-body in- 
teractions (Bl), 4-body interactions (B2) and 6-body in- 
teractions (B3). 

By contrast, the nearest- neighbor approximation (C) 
yields the results which considerably deviates from the re- 
sults of the direct method quantitatively, although some 
qualitative features are still captured. 

In order to investigate the behavior of larger lattices 
L > 35, we have to rely on the approximate methods 
(B1)-(B3) and (C). The domain-wall energy AEp ap and 
the rate of the SW sample r$ w calculated for larger lat- 
tices L > 35 by the methods (B1)-(B3) and (C) are shown 
in Figs. 3(b) and 4(b), respectively. In both methods, 
the ground state under the given BC is searched for be- 
tween the chiral domain-wall state and the 7T-SW state 
by comparing the energies of these two states. As can 
be seen from Fig. 4(b), the rate of the SW sample rsw, 
once increased with L at L < 35, decreases with further 
increasing L and tends to zero in the L — ► oo limit but 
with peculiar oscillations. Likewise, as can be seen from 
Fig. 3(b), AE^ap yields a slope close to 1.9 for large 
enough L, consistent with the value of Ref. [2]]]. Thus, 
the asymptotic large-L behavior seems consistent with 
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FIG. 3: The domain-wall energy A.Ep,ap of even samples, 
calculated by the three methods (A), (B1)-(B3) and (C) men- 
tioned in the text, are plotted versus L on a log- log plot. The 
dotted line is the power law L" 1 - 899 ' • predicted in Ref. jgj. 



Ref. [21| , while the approach to the large-L asymptote is 
rather slow, realized only for lattices with L ^ 40. 

Of course, the methods (B1)-(B3) and (C) assume 
properties of the ground-state configurations are not 
completely rigorous. However, the fact that the non- 
trivial (non-monotonic) small-L behavior revealed by the 
direct method (A) is also reproduced by these approxi- 
mate methods gives some credence to the reliability of 
the approximate methods even for larger L where the 
direct method is not available. 

In Fig. 5, we show the distribution function of the 
domain-wall energy AEj^ap for the sizes of L = 35, Fig. 
5 (a), and L = 240 — 960, Fig. 5 (b). As can be seen 
from Fig. 5(a), the SW samples contributes the com- 
ponent near the edge of the distribution whose weight 
decreases with increasing L. Interestingly, the distribu- 
tion is not smooth at all. As L increases, more bands 
of spikes appear closer to the center AEp^ap = and 
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FIG. 4: The L-dependence of the rate of the spin-wave sam- 
ples within the even samples under the periodic boundary 
condition, calculated by the three methods (A), (B1)-(B3) 
and (C) mentioned in the text. 



FIG. 5: The distribution function of the domain- wall energy 
ABp^p for even samples with L = 35 obtained by the method 
(A) (a), and L = 240,480,960 obtained by the method (Bl) 
(b). 



the amplitude of the spikes at smaller |ASp t Ap| becomes 
larger. Presumably the bands of the spikes reflect nearly- 
discrete spectrum of the distribution of effective inter- 
actions between chiralities, which is not explicit in Eq. 
©■ 

3.1.2 The domain- wall energy: AE C 

Now, we turn to the second type of domain-wall en- 
ergy, AE C , the absolute value of the ground-state energy 
difference between under the R and min(P,AP) BCs. Un- 
der the R BC, the sign of the chirality is reversed at the 
boundary. Therefore, if the sign change in the chiral- 
overlap Oi — mi(R)mi(min(P, AP)) occurs at the bound- 
ary, it actually means that there is no chiral domain- 
wall at the boundary. With this understanding, under 
R/min(P,AP) of even samples, a single chiral domain- 
wall should be introduced into the sample, not accompa- 
nying the SW excitation[2l|. Hence, AE C should be char- 
acterized by the chiral stiffness exponent y K = 1.899 



Indeed, our direct method (A) has fully confirmed this 
expectation. 

In Fig. 6(a), we show the size dependence of the 
domain- wall energy, AE C , estimated by the three meth- 
ods (A), (B1)-(B3) and (C). As mentioned in §3.1.1, in 
the method (B1)-(B3), the position of a chiral domain 
wall is determined by calculating and comparing the en- 
ergies corresponding to all possible Nf T positions of a chi- 
ral domain wall. In the size range L < 35 where the direct 
calculation is available, the data yield a slope about 1.58, 
which is considerably smaller than the expected value 
y K = 1.899 • • •. However, the data for larger lattices ob- 
tained by the approximate methods (B1)-(B3) and (C) 
yield an asymptotic slope consistent with the expected 
value y K = 1.899 • • •. Again, the approach to the asymp- 
totic behavior turns out to be rather slow. 

3.2 Odd samples 
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FIG. 6: The domain-wall energy AE C of even samples, calcu- 
lated by the three methods (A), (B1)-(B3) and (C) mentioned 
in the text, are plotted versus L on a log- log plot. The dotted 
line is the power law ^ _1899 - ■ predicted in Ref. |2l|. 



In this subsection, we deal with the other subset of 
samples where the total number of frustrated plaquettes 
is odd. 

3.2.1 The domain-wall e nerg y: A-Ep^p 

For odd samples, Ref. [21j shows that the P and AP 
BCs always yield exactly the same ground-state energy, 
i.e., AE^ap = 0. In the ground-state chirality pattern, 
a single misfit is introduced into the alternating chirality 
pattern, but always in the same position between under 
the P and AP BCs so that there is no chiral domain-wall 
between under the P and AP BCs. Instead, the 7T-SW is 
generated between them. Indeed, we have confirmed this 
expectation by the direct method (A). 

3.2.2 The domain-wall energy: AE C 

In the case of R/min(P,AP) of odd samples, a single 
chiral domain-wall is expected to be introduced with ac- 
companying the 7r/2-SW[2]]]. This has been confirmed 



by our direct method (A). A typical example of the spin 
configurations under the R and P BCs is shown in Fig. 7. 
The existence of the 7r/2— SW is clearly visible. 



2leg ladder ±J (odd) L=25 



Reflecting 
Periodic 



2leg ladder ±J (odd) L=25 



R-P 



^ ^ * 



FIG. 7: A typical example of the ground-state spin configura- 
tion of odd samples under the periodic and reflecting bound- 
ary conditions, between which a spin- wave of a turn angle n/2 
is generated. The lattice size is L = 25. In the lower panel, 
the relative deviation angle between the spin directions under 
the P and R BCs is illustrated by arrows. 

Here, the domain- wall energy, AE C , is expected to be a 
sum of the chiral domain- wall contribution characterized 
by the chiral exponent y K = 1.899 • • • and the SW con- 
tribution characterized by the SW exponent y' s = l[2lj|. 
For large enough L, a slowly-decaying component, i.e., 
the SW component, should dominate the asymptotic be- 
havior of AE C . In Fig. 8, we show on a log-log plot the L- 
dependence of the domain- wall energy, AE C , calculated 
by our methods (A), (B1)-(B3) and (C). In the methods 
(B1)-(B3), the position of a single chiral domain wall is 
determined by calculating and comparing the energies 
corresponding to all possible N{ T positions of a chiral 
domain wall. In the size range L < 35 where the di- 
rect calculation is available, the data yield a slope about 
1.07, which is slightly larger than the expected asymp- 
totic value y K = 1. This deviation for smaller L might 
be due to the residual contribution of the chiral domain- 
wall. Meanwhile, the data for larger lattices obtained by 
the approximate methods (B1)-(B3) and (C) yield the 
asymptotic slope fully consistent with the expected value, 

y = i. 

3.3 All samples 

As shown above, the behaviors of the domain- wall en- 
ergy largely differ between in even and odd samples. The 
behavior observed when one measures the domain-wall 
energy averaged over all samples can be obtained imme- 
diately by simply combining the above results for even 
and odd samples with equal weights. 
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FIG. 8: The domain-wall energy AE C of odd samples, calcu- 
lated by the three methods (A), (B1)-(B3) and (C) mentioned 
in the text, are plotted versus L on a log- log plot. 



The behavior of AEj^ap averaged over all samples is 
exactly the same as that for even samples, since AEp^p 
is identically zero for odd samples. It exhibits an asymp- 
totic behavior characterized by the stiffness exponent 
y K = 1.899 • • • for large enough L, while there is a signif- 
icant finite-size correction and the asymptotic behavior 
sets in only at L c 40. 

By contrast, AE C is a sum of the chiral contribution 
and the SW contribution, the large-L behavior being 
dominated by the latter. Hence, AE C averaged over all 
samples is asymptotically characterized by the SW expo- 
nent y' s = 1, in contrast to A-Epap- 

These asymptotic behaviors were just as predicted by 
Ney-Nifie et al H2. 



§4. Numerical results of the correlation length 

In this section, we numerically investigate the tem- 
perature dependence of the spin and the chiral corre- 



lation lengths, £ s and £ K , with interest in the associated 
correlation-length exponents, v s and v K . Since the corre- 
lation functions arc bulk quantities, one usually believes 
that they are independent of the type of BCs, or whether 
the sample is either even or odd, in contrast to the case 
of the domain- wall energies analyzed in the previous sec- 
tion. Indeed, this was implicit in the analysis of Ref . |2l| . 

The two-point spin-spin and chirality-chirality correla- 
tion functions, C S (R) and C K {R), are defined by 



C S (R) - [< 5(0,1) • S( RA ) > ], 
C K {R) = [< k kr > 2 ], 



(10) 

(11) 
1 de- 



where < • • • > denotes the thermal average and [• 
notes the average over the bond disorder. 

We directly calculate these correlation functions of the 
cosine model under the P BC by means of the stan- 
dard Monte Carlo simulation. The spin and the chiral 
correlation lengths are extracted by fitting the calcu- 
lated correlation functions by a simple exponential form, 
Aexp[— (r/£)] (A is a constant). To guarantee that the 
estimated correlation lengths are free from the finite-size 
effect, the data are limited to the temperature region 
T/J> 0.1 where both correlation lengths £ s and £ K are 
much smaller than the system size L = 100. At the tem- 
perature T/J = 0.1, £ s and £ K become around 6 and 
3, respectively. In order to be sure that the correlation 
lengths are insensitive to whether the sample is even or 
odd, we estimate £ s and £ K for each case of even and odd 
samples. As expected, £ s of even and odd samples agree 
within the error bars, so does £ K . 

In Fig. 9 (a), we show on a log- log plot the temper- 
ature dependence of £ s and £ K averaged over all sam- 
ples, total number of samples being 100. As can be 
seen from Fig. 9 (a), in the investigated temperature range 
T/J > 0.1, the data yield a slope close to unity for both 
£ s and £ K . This value, unity, may be related to the SW 
stiffness exponent y' s = 1 via the relation v = 1/y. It 
is significantly smaller than the asymptotic value ob- 
tained in Ref.JU v s — v K = l/y K — 0.5263 •• •. How- 
ever, since the temperature range covered in the present 
simulation is rather high and the correlation lengths still 
stayed shorter than the crossover length ~ 40 estimated 
in the previous section, it is quite probable that we need 
to go to lower temperatures to see the true asymptotic 
critical behavior associated with the T — transition. 
Unfortunately, we cannot directly evaluate the correla- 
tion lengths in this low temperature region because of 
the finite-size effect and the thermalization problem. 

Ney-Nifie et al gave an analytic expression of the spin 
and the chiral correlation lengths in the thermodynamic 
limit on the basis of the two assumptions mentioned 
above |2l). The chiral correlation function is obtained 
as 



C K (R) = [<2m -2m R > 2 } 



(12) 
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FIG. 9: The temperature dependence of the spin and the 
chiral correlation lengths, £ s and £ K , of the bulk sample, esti- 
mated by Monte Carlo simulation and numerical analysis of 
the effective charge model. 



s 

J| tanh 2 

i=i 
-R/(* 



k B T 



(13) 
(14) 



where T is the temperature. In Eq. the product is 

taken over all nearest-neighbor effective bonds Vj which 
lie in the section between the sites and R and are la- 
beled as j = 1,2, ... ,S where S is the number of frus- 
trated plaquettes between the sites and R. The chiral 
correlation length £ K behaves in the T — ► limit as 



l 



l/y K = 0.5263. 



(15) 



In the Villain model, the spin correlation function is ob- 
tained as 

C S (R) = [< cos(Vi) - Vi))> 2 ] ( 16 ) 
= e- R ^C chm . gc (R), (17) 



where the SW correlation length £' s is given by 

e = i 



(18) 



Here the spin correlation function is factorized into the 
two parts, one due to the SW and the other due to 
the charges C c harge- Within the effective model studied 
by Ney-Nifle et al [§lj , i.e. the Ising model with the 
nearest-neighbor random antiferromagnetic interaction, 
the charge part becomes, 



Ccharge \R^) 



5/2 / V \ 



is=i 



(19) 
(20) 



In the cases where S is an odd integer, the charge (chi- 
ral) correlation function vanishes. These results do not 
depend on the type of BCs nor on whether the sample is 
either even or odd. 

In Fig. 9, we show these analytic results of the spin and 
chiral correlation lengths for the effective model together 
with the corresponding MC estimates for the original XY 
SG model. In our analytic calculations, we evaluated the 
averages over the disorder in Ea. H13|) and Ea. H19(l tak- 
ing into account the true discrete spectrum of the distri- 
bution of spacings between frustrated plaquettes rather 
than using the continuous expression Eq.Q. As can 
be seen from Fig 9(a), in the higher temperature range 
T/J > 0.1, the analytic results of Ref.|2lJ agree with 
our MC results, exhibiting the near 1/T-behavior. Such 
an agreement observed at higher temperatures might give 
some credence to the reliability of the approximate meth- 
ods. At lower temperatures where the MC result is no 
longer available, the analytic results of Ref.|2l| tend to 
level off, exhibiting a clear crossover. There, for both 
cases of £ s and £ K , a power-law behavior with a much 
smaller asymptotic exponent v s = v K = 0.5263 • • • are 
eventually realized. The crossover from the 1/T behav- 
ior to the i/T 5263 -- behavior occurs below T/J ~ 0.1, at 
the length scale of L = 30 lattice spacings. This crossover 
might be related to the domain-wall result in the previous 
section where a crossover takes place at around L = 40. 

Hence, the asymptotic critical behavior of the spin and 
chiral correlation lengths sets in only at low tempera- 
tures T/J 5 0.1 and at longer length scale £ ^ 30. At 
higher temperatures T/J ^ 0.1 and at shorter length 
scale £ IS 30, a different power-law behavior with an ap- 
parent exponent v ~ 1 fits the data better. 

It might be worth emphasizing here again that, al- 
though the spin and the chiral correlation lengths ex- 
hibit the same critical behavior, apparently with only 
one diverging length scale, there in fact exist two dis- 
tinct length scales at the T = transition of this model. 
This has been already evident in Eqs. (17) and (20), 
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where the spin correlation function is written as a prod- 
uct of the Z2 chiral part, characterized by the correlation 
length with the chiral exponent v K = 0.5263 • • • and the 
50(2) SW part, characterized by the correlation length 
with the SW exponent v' s = 1. Hence, while there actu- 
ally exist the two diverging length scales at the T = 
transition of the model, the one diverging more slowly, 
i.e., the one with smaller v, dominate the asymptotic be- 
havior of spin correlations, masking the existence of the 
other correlation length which diverges more rapidly. 

We note that such a "masking" phenomenon arises 
only when the inequality v K < v s holds between the Zi 
and 50(2) correlation- length exponents. If this inequal- 
ity would be opposite, the masking phenomenon would 
not show up in the spin correlations. Then, the existence 
of two correlation lengths would manifest itself more di- 
rectly in the associated correlation functions, the 50(2) 
spin correlation length in the spin correlations and the 
Zi chiral correlation length in the chiral correlations. 
This actually occurs in the aforementioned regularly frus- 
trated ID XY model, where one has v K = 00 > v s = 1 
as shown rigorously for the case of ID triangular lattice 
inRef. 0. 

Naturally, one expects essentially the same behavior in 
the present two-leg ladder XY model. Let us consider the 
regularly frustrated two-leg ladder XY model such that 
all spacings between frustrated plaquettes are equal to I. 
The analytic expressions of the charge (chirality) and the 
spin correlation functions, Eqs.(13) and (19), expected to 
be valid at low enough temperatures T/J <C 1, can also 
be used in the regular case where the average over the 
bond disorder should be dropped. As examples, we show 
in Fig. 10 the temperature dependence of the spin and the 
chiral correlation lengths £ s and £ K for the cases of I = 2 
and 4. As can be seen from the figure, the chiral corre- 
lation length £ K outgrows the spin correlation length £ s 
at some finite temperature T x . At higher temperatures 
T > T Xl £ s and £ K exhibit a more or less similar behav- 
ior with £ s > £ K . At lower temperatures T < T x , while 
£ K continues to exhibit a behavior similar to the behav- 
ior observed at higher temperatures T > T x , the spin 
correlation length £ s dramatically changes its behavior. 
The growth of £ s with the decrease of the temperature is 
dramatically slowed down below T x yielding £ s <C £ K ; a 
manifestation of the spin-chirality decoupling. 



§5. Summary and discussion 

We numerically investigated the domain-wall energies 
and the spin and the chiral correlation lengths of the ID 
±J XY SG ladder. Analytic results obtained by Ney- 
Nifle et al were confirmed for asymptotically large L, 
while the finite-size effect could be significant. Concern- 
ing the domain-wall energies, the asymptotic behavior 
sets in only for lattices with L £ 40. Concerning the 
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FIG. 10: The temperature dependence of the spin and the 
chiral correlation lengths, £ s and £ K , of the two- leg ladder reg- 
ularly frustrated XY model evaluated by the effective charge 
model. The spacings between frustrated plaquettes are equal 
to /. Here the results of I — 2 and I — 4 are shown as exam- 
ples. 

correlation lengths, the asymptotic behavior sets in only 
at low temperatures T/J % 0.1. 

The domain-wall energies show different behaviors de- 
pending on the type of BCs and whether the number 
of frustrated plaquettes is even or odd. The domain- 
wall energy associated with the P/AP BCs, ALp.ap, is 
identically zero for odd samples, and exhibits a rather 
complex behavior for even samples. For even samples, 
although the lowest-energy excitation is always a chiral- 
domain wall pair for asymptotically large L, the one for 
smaller L is sometimes a 7T-SW without accompanying 
the chiral domain wall, which eventually gives way to 
the chiral domain- wall pair excitation for larger L. Some- 
what unexpectedly, the approach to the asymptotic large- 
L limit is rather slow and could be even non-monotonic 
(see Fig. 4(b)), and the asymptotic large-L behavior sets 
in only at L c 40. The chiral domain wall and the SW 
bear the associated stiffness exponents, y K — 1.899- •• 
and y' s = 1, respectively. We note that, if the inequality 
between y K and y' s would be opposite, i.e., if y K < y' s , 
then A£p,ap would be characterized by y s , not by y K . 

The domain-wall energy associated with the 
R/min(P,AP) BC, AE C , is governed by a chiral 
domain-wall excitation, with and without the SW for 
odd and even samples. Absence or presence of the 
SW excitation gives rise to the different asymptotic 
behaviors of AE C for even and odd samples, respectively. 
The SW excitation, whenever it is induced, governs the 
large-L asymptotic behavior of the domain wall energy 
since y' s = 1 < y K = 1.899- •• in the present model. 
Thus, AE C is characterized by an asymptotic stiffness 
exponent y — y K and y — 1 for even and odd samples, 
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respectively. We emphasize again that, if the inequality 
between y K and y' s would be opposite, AE C would always 
be characterized by the chiral exponent y K for both even 
and odd samples. 

We also numerically investigated the behavior of the 
spin and the chiral correlation lengths. Both exhibit 
the divergence characterized by the chiral exponent v K = 
1/z/k = 0.5263 • • •, whereas this asymptotic behavior sets 
in only at low temperatures T/J % 0.1 and at longer 
length scale £ > 30. At higher temperatures T/J > 0.1 
and at shorter length scale £ |$ 30, a different power-law 
behavior characterized by an apparent exponent v ~ 1 
is realized. Although both the spin and the chiral corre- 
lation lengths exhibit the critical behavior with a com- 
mon exponent y K , the system in fact possesses two dis- 
tinct diverging length scales: one associated with the 
50(2) SW and the other associated with the Z 2 chi- 
rality. In that sense, the SO(2) part and the Zi part 
are decoupled in this model. Reflecting the inequality 
v' a = l/y' s = 1 > v K = l/y K = 0.5263- • •, however, the 
spin correlation is dominated by the chiral exponent v K , 
not by the SW exponent v' s = 1. The inherent spin- 
chirality decoupling of the present model is then masked. 

These behaviors of the domain- wall energies and of the 
correlation lengths of the ID ± J XY SG ladder are seem- 
ingly at odds with the behaviors of the corresponding 2D 
and 3D models suggested in Refs 0, 0|. In the 2D 
and 3D XY SG, Refs. ([Hj) conjectured that AE P)AP 
was governed by the SO(2) spin exponent, while AE C 
was governed by the Zi chiral exponent. This should be 
contrasted with the behavior of the present ID model 
where AiSp^AP was governed by the Zi chiral exponent 
while AE C was governed by the SO(2) spin (SW) expo- 
nent. As argued in the previous section, however, this 
apparent difference is simply the consequence of the dif- 
ference in the relative magnitude of the Z2 and 5*0(2) 
exponents between the two models. Indeed, in the 2D 
case, Ref . [l(ij| estimated y K ~ 0.5 and y s ~ 1.0 from the 
L-dependence of AE C and AEp^p, respectively, where 
one has y K < y s in contrast to the present ID case. If 
the inequality between the two stiffness exponents y K and 
y' s were opposite in the present ID model, i.e., y K < y' s , 
essentially the same behavior as that in the 2D and 3D 
models, i.e., AE C L~ v ° and AEp^xp ~ L~ Vk , would 
arise. Hence, the apparent large deviation from Ref. [13 
is merely due to the difference in the relative magnitudes 
of the two stiffness (and correlation-length) exponents. 

In this connection, a significant difference between 
the present ID ladder model and the higher-dimensional 
models might be that, while the effective interaction be- 
tween the charge variables m, is short-ranged in ID, that 
in higher dimensions is long-ranged. This assigns a non- 
trivial character to the charged excitation in the higher 
dimensional models. In 2D, this charged excitation could 
be viewed as a vortex. (In 3D, it is a vortex-line.) The im- 



portant characteristic of such a vortex excitation is that it 
breaks the charge neutrality condition and interacts with 
each other via the long-ranged Coulombic interaction. In 
the present ID ladder model, even though there certainly 
exists an excitation which apparently breaks the charge- 
neutrality condition, the associated interaction between 
them is always short-ranged, and cannot be regarded as 
a genuinely-charged vortex-like excitation. 

Thus, in higher dimensions, excitations associated with 
the charge excitation are at least of two distinct types: 
The one is the vortex excitation which breaks the charge- 
neutrality, and the second one is the ordinary chiral ex- 
citation which preserve the charge- neutrality. In 2D, the 
former one is a point defect (vortex), while the latter 
is a line defect (domain- wall) . Note that the former 
vortex-like excitation exists in both frustrated and un- 
frustrated XY spin models, often playing a vital role in 
the order-disorder process, while the latter chiral excita- 
tion is peculiar to the frustrated XY spin models. We 
expect that, in higher dimensions, this genuinely-charged 
vortex excitation might give rise to the another diverging 
length scale which predominantly disorder the spin (not 
the chiral) order. A possible conjecture would be that the 
stiffness exponent associated with the vortex excitation 
y v takes a value smaller than the SW stiffness exponent 
or the chiral stiffness exponent, and is more effective in 
disordering the spin. If this is really the case, the spin 
correlation length and the domain-wall energy AEp^p 
would be governed by the non-chiral exponent y v asso- 
ciated with the genuinely-charged excitation. In the 2D 
XY SG with the ± J interaction, the numerical estimate 
of Refs.[10| suggested y v ~ 1, which appears to be larger 
than y K ~ 0.5. 

We note that, even in ID, the long-ranged interaction 
between the charge variables could arise in some special 
cases, e.g., in the tube lattice investigated by Hill et al 
|22|. In the tube lattice, however, the application of nei- 
ther the P, AP nor the R BCs is capable of generating the 
genuinely-charged vortex excitation, i.e., the excitation 
in the q + variable in the notation of Ref.|22|. In fact, the 
application of either the P, AP and R BCs generates only 
the charge-neutral chiral excitation, the excitation in the 
q~ variable in the notation of Ref. [2^], which interacts 
only via the short-ranged interaction. As such, the role 
of the genuinely-charged vortex excitation still remains 
to be seen in the ID tube model. In order to elucidate 
the nature of the spin and the chirality orderings of the 
frustrated XY"spin systems, it would be important to fur- 
ther clarify the role of the vortex-like genuinely-charged 
excitation, not only in the ID tube model but also in the 
higher-dimensional models. 

In the present paper, we concentrated on the XY 
SG ladder with the ±J interaction. The corresponding 
model with the Gaussian interaction was also studied in 
the literature, e.g., in Ref. |23| . Our numerical study sug- 
gested that these two models, i.e., the ID XY"SG ladder 
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with either the ±J or the Gaussian interaction exhibit 
quite different behaviors. The properties of the ID XY 
SG ladder with the Gaussian interaction will be reported 
elsewhere !24l. 



APPENDIX: EFFECTIVE CHARGE 
HAMILTONIAN 

In this appendix, we briefly summarize the deriva- 
tion of the effective charge Hamiltonian corresponding to 
the original XY spin-glass model described by the spin 
Hamiltonian Eq. QJ. We note that the Villain's Hamil- 
tonian, which contains the two-body interaction between 
the charges only, is not exact even in the T — * limit, 
contrary to a common belief. In the following, we give 
explicit forms of the correction terms to the Villain's ap- 
proximation, which includes certain 4, 6,... body effective 
interactions between the charges. Detailes will be re- 
ported elsewhere wtihin a more general context |2fij . 

As usual, the starting point is the identity 



= /3 cos 6 



E e*%G8), 



(21) 



by the two external change variables m ex i and m 0X 2 no- 
ticed in [2l]]. In the case of periodic (P) and antiperi- 
odic (AP) boundary conditions (BC), all charge vari- 
ables are subjected to the global neutrality condition 
J2i m i + ^Bxj + ?7iox2 = due to the presence of a mass- 
less mode |2l|. In the notaion of §2, the external charges 
reads as m cx i = n — ^ for P BC and m cx i = n — 
for AP BC while m GX 2 = rrij — m cx i. One finds that 
the same global neutrality condition applies also for the 
higher-body terms given below for the cases of P and AP 
BCs. In the case of reflecting (R) boundary condition, on 
the other hand, one finds to ox i = to CX 2 = in the dual 
mapping and also finds that the massless mode and thus 
the global neurtarily condition is absent. 

The explicit form of the higher-body terms becomes 
increasingly complicated. The 4-body Hamitonian reads 
as, 



/4— body 



W 4 - body = E [Df + (D exl )j + (D cx2 )t] , (25) 



with 



u, 



u, 



(26) 



where j3 — 1/T is the inverse temperature and I p (f3) is 
a modifeid Bessel function. Now we use a useful formula 
for the assymptotic behaviour in the limit (3 



00, 



log W)/- 



\ - c, 



<«■ rp. 



2m-l 



(l + 0(T))/i" 



(22) 

where /i = Ap 2 . Here terms which vanish in the T — > 
are repsented as 0(T). The first few coefficients reads as 
d = -1/8, c 2 = 1/192, c 3 = -3/2560, c 4 = 15/28672, 
c 5 = -35/98304, c 6 = 945/2883584,... Note that the 
usual Villain's approximation amounts to assume that 
C\ = —1/8 and c m = for m > 1. 

Following the standard steps of mapping of the ori gina l 
spin model to the charge model on the dual lattice [25j, 
one finds that the effective charge Hamiltonian in the 
T — > limit can be written as, 

^charge = H2-body + ^4-body + H6-body--. (23) 

The terms on the r.h.s are due to the terms of m = 
1, 2, 3, ... in Eqs. The first term is nothing but the 

usual Villain's Hamiltonian which describes the two-body 
charge interactions, 

^2-body = £ U^mj + ^-(m c 2 xl + ml,), (24) 



and 



(-De X l)i 



(A, 



2tt 



m ° xi + E 



-nij, 



2vr 



L 



E 



Ui 



The 6-body Hamiltonian reads as 
4 3 c 3 



(27) 



(28) 



•^6 — body 



3! 



4 m i M 



- E(^j + B ij) D i (Pexl)f + (D ex2 )!) 



i 

L 



£(a, 



E( Dex2 )i 



where 



_ 1 cos(fc(i - j))(l - cos(fc)) 



2 — cos fe 



(29) 



(30) 



where [7jj is given by Eq. JSJ. The second term on 
the r.h.s describes global spin-wave excitations induced 



B: 



1 ^ r sin(fc(i — j)) sin(fc) 



E 



2 — cos k 



(31) 
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_ 1 ^ cos(fc(i 
ij ~ ~~L 4-*' 4 - 2 cos fc 



i)) 



(32) 



For a simple demonstration, let us consider a special 
2-leg ladder sample in which all plaquettes are unfrus- 
trated. In this particular sample, one knows that a spin- 
wave is induced under the AP boundary condition so that 



AE 



PAP 



= -2L 





- 1 


cos (l) 







7T~ 

T 



21 L 3 



1 7T 6 

6! Is 



(33) 

exactly. One can check that the first three terms in the 
last equation can be obtained by the energies associated 
with the 2-body (Villain's approximation), the 4-body 
and the 6-body interactions, respectivly. 
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